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Abstract 

In y^iivjv = 200 GeV Au+Au collisions PHENIX reported a significant enhancement in the low-mass region 
(0.1 < m ee < 0.7 GeV) of the dielectron spectrum, which is still not fully understood. Several theoretical 
works and an indirect measurement suggest that, due to the possible restoration of the £/a(1) part of the 
chiral symmetry in a hot and dense medium, the mass of the rf meson may substantially decrease. This 
work reports on a statistically acceptable description of the PHENIX low-mass dilepton enhancement using 
a radial flow dominated meson spectra, chain decays of long-lived resonances and an in-medium rf mass 
modification. 
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1. Motivation 

In comparison to the expected hadronic decays, 
a significant excess has been measured in the (0.1 — 
0.7) GeV invariant electron-positron mass range of 
the PHENIX dilepton cocktail 0. To quantify 
this deviation, the corresponding \ 2 an d confidence 
level (CL) were calculated; x 2 /NDF = 50.3/12, CL 
= 0.01% (for 0.12-1.2 GeV mass range). Prelim- 
inary results from STAR also show a substantial, 
although smaller effect in a different geometrical 
acceptance . Although this low-mass dilepton en- 
hancement is a subject of detailed theoretical inves- 
tigations, according to a recent conference report 0] 
this enhancement is not yet fully understood. 

Lattice results Q suggest that the transition 
to chirally symmetric and deconfined matter is a 
smooth, cross-over type transition, where the chiral 
symmetry restoration as measured by chiral suscep- 
tibility happens at a significantly lower temperature 
than the confinement of quarks to hadrons. This 
was measured by the Polyakov loop: the continuum 
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extrapolated value for the transition temperature of 
the renormalized Polyakov loop is 25+4 MeV higher 
than that of the chiral susceptibility. This result 
suggests an interesting possibility of chiral hadron 
dynamics in the 150-175 MeV temperature range. 
Chiral symmetry restoration may change various 
properties of hadrons, such as their mass, width, 
branching ratios and decays. 

The PHENIX collaboration was the first to re- 
port on a measurement of direct photon spectrum 
in y/SNN = 200 GeV Au+Au collisions in Ref. @, 
where the average temperature of the photon emit- 
ted region was found to be 219+19 (stat)±19 (syst) 
MeV, implying that the initial temperature was at 
least about 300 MeV, and the (kinetic) freeze-out 
temperature of the order of 140-170 MeV. 

The in-medium drop of the mass of the rf me- 
son was anticipated long time ago @, 0, B- The 
rf meson is different from the other pseudoscalar 
mesons in the sense that due to the axial anomaly it 
has a very large mass (m n i = 958 MeV) which pre- 
vents it from acting as the 9 th pseudo-Goldstone bo- 
son. However, in a hot and dense hadronic medium 
where the chiral symmetry is partially restored, its 
mass can be reduced to near its quark model value. 

The first indirect observation of in-medium chi- 
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ral dynamics, based on Bosc-Einstcin correlation 
(BEC) measurements and on the analysis of rf — > 
ir + ir~7] 7r + 7r _ (7r + 7r _ 7r°) decay chain, suggested a 
significant and an at least 200 MeV in-medium mass 
drop of the 77' mesons in these reactions [i|,[Kl 11 1. 
These results triggered renewed interest in theo- 
retical investigations of Ua(1) symmetry restora- 
tion: the Witten-Veneziano formula has been gen- 
eralized to finite temperature to check if the in- 
directly measured in-medium rf masses are in the 



theoretically acceptable range [12j, lattice QCD cal- 
culations were performed to study the temperature 
dependence of the Ua(^) symmetry restoration [l3| . 
nuclear bound states and partial restoration of the 
77' mass was investigated in Ref. 14|, just to men- 
tion the some of the most inspiring new results. In 
this manuscript we cross-check this indirect obser- 
vation, based on a hadronic decay channel, utiliz- 
ing the dilepton (Dalitz) decay channel of the rj 
mesons. 

We have to keep in mind that changes in other 
meson properties may also occur in the hot and 
dense hadronic medium. For instance the broaden- 
ing of the Vj width, which was reported in a recent 
measurement of the TAPS collaboration [l5[ in cold 
nuclear matter. Nuclear photoproduction calcula- 
tions on the p meson suggest a subtantial broaden- 
ing which may account for part of the excess low- 
mass dileptons [lf|. The NA60 Collaboration re- 
ported the first measurement of a p spectral func- 
tion, based on a low-mass dilepton measurement 
in a similar kinematic range as the PHENIX low 
mass dilepton enhancement, in £x a b = 158 AGeV 
indium-indium collisions at the CERN SPS [17|. 
The associated spectral function of the p meson is 
found to indicate significant broadening but within 
experimental errors vanishing in-medium p mass 
modification [l7| . The same mechanism, however, 
is not able to fully account for the PHENIX low- 
mass dilepton enhancement in y/sJrN = 200 GeV 
Au+Au reactions. Also note that, contrary to 
RHIC collisions, a previous analysis of NA44 S+Pb 
data in the BEC channel showed no significant 77' 
mass drop at CERN SPS energies [11]. This indi- 
cates that the chiral dynamics is different at RHIC 
energies of ^/saT/v = 60, 130 and 200 GeV from 
the observations at the CERN SPS top energy of 
y/spjN = 19.4 GeV (see Ref. [6] for details on this 
point). 

The method we utilize here is fully consistent 
with the methods we applied in Refs. [7,8] to ob- 
tain the indirect 77' mass modification measurement 



using Bosc-Einstcin correlations of pions from 7/ 
chain decays. As far as we know (also see Ref. Q), 
this is the first statistically acceptable description 
of the low-mass enhancement in 200 GeV Au+Au 
collisions at RHIC, as we shall demonstrate below. 
However there are several possibilities for further 
improvements. Ideally all the in-medium effects 
should be considered simultaneously. In this anal- 
ysis, however, we restrict ourselves in exploring the 
complex phenomenon of chiral symmetry restora- 
tion only to the mass modification of the 77' meson, 
which is expected to be the most prominent effect 
especially at the lower temperatures of the hadronic 
phase. 

The Dalitz-dccay (7/ — > e + e~7) of the 77' meson 
contributes to the low mass region (below ~1 GeV), 
where the observed excess is present [l| . 

2. The rf meson in a hot, dense medium 

2.1. Effect of mass-modification 

As a result of a possible chiral symmetry restora- 
tion the mass of the 77' meson can significantly 
change, but this mass modification can not be ob- 
served directly. Since the 77' meson decays in the 
vacuum with its original mass and it has a rel- 
atively small interaction cross-section, its lifetime 
is long, compared to the fireball. Thus we can 
search for this effect via its enhanced production, 
which can be described with a modified Hagcdorn- 
formula QfTl: 



Srf = 



(1) 



where m v i is the vacuum mass- 



m*i is the in- 



medium mass of the 77' meson, while T conc i is the 
characteristic slope of the 77' mesons within the 
formed medium. In order to explore the role of 
77' meson enhancement on the -^/SaT/v = 200 GeV 
Au+Au collision data of PHENIX 0, extended 
EXODUS simulations were elaborated to generate 
different dilepton spectra corresponding to physi- 
cally reasonable 77' properties. 

2.2. The shape of the 77' spectrum 

According to the conservation of energy, those 
77' mesons that can escape from the hot and dense 
hadronic medium, while regaining their mass, must 
lose momentum. Denoting the properties of the in- 
medium 77' with an asterisk, while quantities with- 
out further notation correspond to the free 77', the 
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Figure 1: The reproduced PHENIX dilepton spec- 
trum of Ref. P and a resonance pr < 2 GcV cut. 
The bottom shows the ratio of the two, proving the 
dominance of the low-px cocktail. 



conservation of energy requires rf mesons to fulfill 
the 



"■V rP T =™-,+Pt (2) 

relation. This condition results in two different 
components of the rf spectrum. Some of the rf 
meson's transverse momenta is large enough to let 
them escape from the hot and dense medium to the 
asymptotic vacuum by decreasing their momenta 
and increasing their mass. Those rf mesons, how- 
ever, that have low transverse momentun|f| in the 
chirally modified vacuum, cannot come on mass 
shell thus they are trapped in that modified vac- 
uum. They are released and come on-shell only af- 
ter the dissociation of the modified vacuum state to 
the asymptotic vacuum. Both components realize 
an effective |0)* + \r]'}* —> |0) + transition, but 
in a different way. In the first scenario, the trapped 
rf 
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distribution, where B^ 1 is the slope parameter of 
the rf condensate. The high px part of the rf spec- 
trum, however, follows a spectrum similar to other 
hadrons, and can be described with the same hy- 
drodynamical models. 

The final rf spectrum is the sum of Eq. [3] and 
the hydrodynamical spectrum, which results in a 
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double exponential spectrum as used in Refs. [HGiJ 
and also in the current analysis. 

3. Radial flow for dilepton enhancement 

Another idea which might bring us closer to un- 
derstand the excess is to introduce the mT-scaling 
with radial flow of the mesons, which was not in- 
cluded in the orig inal PHENIX analysis Q. Radial 
flow is valid in the low pt range, where it describes 
the particles' spectra well, and this is also the region 
where the dilepton spectrum has the most contribu- 
tion from. On Fig.Q]thc original PHENIX result Q 
and a resonance (pr < 2 GeV) cut result are plot- 
ted, confirming the dominance of the low-p<r part 
of the spectra. Thus examining the effect of the 
radial flow seems to be inevitable, as it might be 
responsible for certain parts of the excess. 

As previous examinations suggested, a hydro- 
motivated spectrum might give a better estimate 
on the low-pT regime for all mesons [Isj . The fol- 
lowing simplified hydro spectrum was chosen: 
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Hydro models usually describe the primordially 
created particles, but not their decay products. 
This effect is most important in case of the ir 
meson, which can be a daughter of many other 
hadrons. Thus we corrected the measured pion 
spectrum using the Core-Halo (C-H) picture [Til 
20l . 21, 22 1 . Primordial pions (together with decay 
products of rapidly decaying resonances that yield 
thermalized pions) are in the core, the inner set of a 
formed medium, while decay product pions of long 
lived resonances form the halo, farther away from 
the collision point. Thus if we divide the measured 
7T° spectrum with 
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we get an estimation for the number of core pions. 
An upper boundary on the corrections for neutral 
pions was obtained by averaging the 7r + and 7r~ C- 
H corrections X H = (A 7r+ + A 7r _)/2, while lower 



boundary was Al 



based on the 



consideration that long lived resonances do not pro- 
duce more than twice as much charged as neutral 
pions. The results are shown in Fig. [21 From now 
on, pion spectra are always C-H corrected. 

The hydro formula of Eq. 0] was simultaneously 
fitted to available ^ , , p and p minimum bias 
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Figure 2: Core-Halo correction parameter, measur- 
ing the ratio of particles in the core compared to the 
total number of particles (including halo). The up- 
per and the lower boundaries of the ir° correction 
band are the average of charged pion corrections 
and half of the average charged pion corrections re- 
spectively. 



data of Ref. [24j, with the same constraints, e.g. the 
fit ranges were 0.2-1.0 GeV/c 2 for tt* and 0.1-1.0 
GeV/c 2 for K ± 1 p and p in (tot — tn). It is shown 
on Fig. |31 as a function of px- 

The fit results are: (u T ) 2 = 0.21 ±0.01, and T Q = 
178 ± 4 MeV, while the a parameter was within 
errors consistent with 0, hence it was fixed to 0. 
The fit quality is: x 2 /NDF = 63/49, CL = 8.6%. 

4. Dilepton cocktail simulations 

4-1. Radial flow and chain decays 

Every meson was generated with the hydro spec- 
trum of Eq. 21 except for the rf. Their normal- 
ization factors, however, could not be determined 
in every case, due to the lack of data of the hydro 
validity range (jpr < 2 GcV). Moreover, these in- 
clusive spectra contain the decay products as well 
besides the core particles. Therefore thermal model 
predictions of Ref. [25[ were used in this work. This 
is a central reasoning of this analysis: from the 
number of core particles we were able to estimate 
the enhanced production of rf mesons, and we can 
conclude the in-medium mass of rf from this en- 
hancement by substituting it into Eq. [TJ 

This different approach requires the introduction 
of chain decays: all the relevant decay chains of 
the mesons contributing to the low-mass dilepton 
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Figure 3: Simultaneous fit of Eq. HJ for n ± K ± , p 
and p data with the constraint from Ref. [24j, e.g. 
the fit ranges are 0.2-1.0 GeV/c 2 for 7r ± and 0.1- 
1.0 GeV/c 2 for K ± 1 p and p in {rax — m). Pions 
are Core-Halo corrected. 



spectrum were taken into account. The simulation 
tool was EXODUS, a Monte-Carlo-based event and 
decay generator already fine-tuned to the PHENIX 
acceptance. 

The simulation results are shown with dashed 
lines on Fig. @] The agreement with data has been 
already improved, the simulated dilepton spectrum 
deviates less from data, the x 2 /NDF = 24.1/12, the 
confidence level CL = 1.9%. 
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Figure 4: Dilepton spectra with the hydro spectrum 
of Eq. |H including chain decays. The dashed lines 
correspond to an original rj' mass scenario (e.g. to 
the hydro background) while the solid lines repre- 
sent the best reduced mass, m*, = 520 MeV. 
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4-. 2. Including rf mass modification 

In the previous subsection, the new, hydrody- 
namical background was evaluated. In this one it 
will be combined with the exotic effect of the 77' 
mass modification. The 77' spectrum is practically 
determined by two parameters: the inverse slope 
parameter, _B _1 (sec Eq. [3]) and the in- medium 
77' mass, 771*/. We scanned through the physically 
meaningful range of these two parameters in order 
to find the most likely setup. A simplification was 
introduced on this step to reduce the simulation 
time: we handled to*, as a parameter which gov- 
erns only the overall size of the 77' excess (see Eq. [T]) , 
while £? -1 is responsible for the shape of the spec- 
trum. In fact, the slope of the longer component 
of the excess is also influenced by the in-mcdium 
mass, but the effect was found to be small in the 
physically interesting in-mcdium mass range. Fig. 
[5] shows the confidence level map of these two pa- 
rameters calculated for the 0.12-1.2 GeV region. 
The confidence level map in [5] indicates the first re- 
sults of the scan of the in-medium 77' mass and the 
slope parameter B~ x . This map is not conclusive 
for the inverse slope parameter B^ 1 in the sense 
that a more subtle analysis, with smaller B~ x bins 
and to*, spectrum-shape dependence is needed to 
determine an accurate value with proper system- 
atic errors. From this preliminary scan, B^ 1 = 55 
MeV can be obtained. In this B~ Y bin a reasonable 
mass range, 480 < to*, < 600 MeV corresponds to 
CL > 15%. Regardless of B~ x , within the observed 
range, to*, < 400 MeV can be ruled out with 99.9% 
confidence level. 

Fig. 0] also shows the best cocktail (B^ 1 = 55 
MeV, m*, = 520 MeV) with a solid line. The agree- 
ment is further improved to x 2 /NDF = 15.2/12, CL 
= 23%. Fig. [6] shows the dilepton cocktail for dif- 
ferent m*, values. 
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5. Conclusions 

The current analyis is aimed at the understand- 
ing of the low-mass enhancement of the PHENIX 
dilepton spectrum of Ref. [l| . A more detailed de- 
scription can be found in Ref. [26[. As a first step, a 
different hadronic cocktail was evaluated based on 
the fact that \ow-px resonances dominate the dilep- 
ton spectrum, as it was examined and its absolute 
normalization detailed in Section [3J 

As far as we know (see Ref. Q ) , this report is the 
first that is able to give a statistically acceptable de- 
scription of the low-mass dilepton enhancement in 
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Figure 5: Preliminary confidence level map of the 
two 77' spectrum parameters; B~ x and 777*,. CL 
> 0.01% is plotted. 



m* 


= 200 MeV 


m* 


= 300 MeV 


m* 


= 400 MeV 


m* 


= 500 MeV 


m* 


= 600 MeV 


m* 


= 700 MeV 


m* 


= 800 MeV 


m* 


= 958 MeV 


200 GeV AuAu 




Figure 6: Dilepton cocktail for different in-mcdium 
77' mass values. 



200 GeV Au+Au collisions at RHIC. Apparently, 
radial flow and chain decays have to be taken into 
account given that this low-mass dilepton enhance- 
ment is in the domain of soft physics with px < 2 
GeV. These effects together provide a data descrip- 
tion with CL=1.9%. Adding in-medium 77' mass 
modification increases the quality of agreement, the 
best values correspond to in-medium 77' mass of 520 
MeV, with a CL = 23%. Extremely low 77' masses 
below 400 MeV can be excluded with 99.9 % confi- 
dence level. The current results summarize the sta- 
tus of our research during August 2012 [27}. In this 
work, the error analysis and other systematic stud- 
ies are still preliminary, and further work is needed 
to investigate in more detail the transverse momen- 
tum and centrality dependence of our findings. 
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